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We develop the theory of discrete time Lagrangian mechanics on Lie groups, originated 
in the work of Veselov and Moser, and the theory of Lagrangian reduction in the discrete 
time setting. The results thus obtained are applied to the investigation of an integrable time 
discretization of a famous integrable system of classical mechanics, - the Lagrange top. We 
recall the derivation of the Euler-Poinsot equations of motion both in the frame moving with 
the body and in the rest frame (the latter ones being less widely known). We find a discrete 
time Lagrange function turning into the known continuous time Lagrangian in the continuous 
• »h ' limit, and elaborate both descriptions of the resulting discrete time system, namely in the body 

frame and in the rest frame. This system naturally inherits Poisson properties of the continuous 
time system, the integrals of motion being deformed. The discrete time Lax representations are 
5/5 , also found. Kirchhoff 's kinetic analogy between elastic curves and motions of the Lagrange top 

is also generalised to the discrete context. 
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1 Introduction 



This paper is devoted to the time discretization of a famous integrable system of classical mechanics 
- the Lagrange top. This is a special case of rotation of a rigid body around a fixed point in 
a homogenious gravitational field, characterized by the following conditions: the rigid body is 
rotationally symmetric, i.e. two of its three principal moments of inertia coincide, and the fixed 
point lies on the axis of rotational symmetry. We present a discretization preserving the integrability 
property, and discuss its rich mechanical and geometrical structure. Notice that until recently |b| 
only the integrable Euler case of the rigid body motion was discretized preserving integrability tyj , 
fljyTvTl , HBLSH - Consult also [|ALl , ||, flDJMfl , |QNCVl for some fundamental early papers on the 
subject of integrable discretizations, and [BP], Q] for reviews of this topic reflecting the viewpoints 
of the present authors and containing extensive bibliography. 

We found the standard presentation of the Lagrange top in mechanical textbooks insufficient in 
several respects, and therefore chose to write this paper in a pedagogical manner, giving a systematic 
account of the new results along with the well known ones represented in a form suitable for our 
present purposes. The paper is organized as follows. 

The introduction recalls the classical Euler-Poinsot equations for the motion of the spinning top 
in the frame moving with the body. Further, we give less known Euler-Poinsot equations describing 
the Lagrange top in the rest frame (they cannot be directly generalized to a general top case). We 
finish the introduction by announcing a beautiful time discretization of the latter equations. 

In order to derive this discretization systematically, we need some formalism of discrete time 
Lagrangian mechanics on Lie groups. The discrete time Lagrangian mechanics were introduced by 



Veselov Q, [ MV ], see also [ WM ], but the case of Lie groups have certain specific features which, 



in our opinion, were not worked out sufficiently. In particular, there lacks a systematic account 
of the discrete time version of the Lagrangian reduction (which is fairly well understood in the 
continuous time setting, cf. [MS], [ HMRfl ). Also, we think that some technical details in @, [MV], 
[WM] could be amended: in working with variational equations these authors systematically use 
Lagrange multipliers instead of introducing proper notions such as Lie derivatives (specific for Lie 
groups as opposed to general manifolds). Therefore we give a detailed exposition of the discrete 
time Lagrangian mechanics on Lie groups in Sect. 3. In order to underline an absolute parallelism 
of its structure with that of the continuous time Lagrangian mechanics, we included in Sect. 2 also 
a presentation of (a fragment of) the latter, which is, of course, by no means original. 

Sect. 4 is devoted to a Lagrangian derivation of equations of motion of the Lagrange top, both 
in the rest and in the body frames. Finally, in Sect. 5 we do the same work for a discrete time 
Lagrange top. 

It has to be mentioned that the actual motivation for the present development came from 
differential geometry, more precisely, from the theory of elastic curves. A brief account of the 
relations between spinning tops and elastic curves is given in Sect. 6. 

We also give three appendices. Appendix A is for fixing the notations of Lie group theory. 
In the Appendix B we collect the main results of Sect. 2, 3 in the form of an easy-to-use table. 
Finally, Appendix C contains some conventions and simple technical results on a specific Lie group 
we work with, namely on SU(2). It should be remarked here that our experience with various 
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integrable discretizations convinced us that working with this group has many advantages when 
compared to the group SO(3), more traditional in this context. 

A standard form of equations of motion describing rotation of a rigid body around a fixed point 
in a homogeneous gravity field is the following: 

M = M x fi(M) + P x A 

(1.1) 

p = p x n(M) 

Here M = (Mi, M 2 , M^) T £ R 3 is the vector of kinetic momentum of the body, expressed in the 
so called body frame. This frame is firmly attached to the body, its origin is in the body's fixed 
point, and its axes coincide with the principal inertia axes of the body. The inertia tensor of the 
body in this frame is diagonal: 

( Ji \ 
J = J 2 . (1.2) 

V j 3 J 

For the vector Q = £l(M) of the angular velocity we have: 

n = fi(M) = J _1 M = (J^Mi, Jz 1 M 2 , J 3 " 1 M 3 ) T G M 3 . (1.3) 

The vector P = (P\, P 2 , Pz) T € M 3 is the unit vector along the gravity field, with respect to the 
body frame. Finally, A = (Ai, A 2 , A^) T G M 3 is the vector pointing from the fixed point to the 
center of mass of the body. It is a constant vector in the body frame. 

It is well known that the system (|l.l| ) is Hamiltonian with respect to the Lie-Poisson bracket 
of the Lie algebra e(3) of the Lie group E(3) of euclidean motions of R 3 , i.e. with respect to the 
Poisson bracket 

{M i ,M J } = e ijk M k , {M i ,P j } = s ijk P k , {P i ,P j } = 0, (1.4) 

where Eij k is the sign of the permutation (ijk) of (123). The Hamiltonian equations of motion for 
an arbitrary Hamilton function H = H(M, P) in the bracket ( |1.4[) read: 

M = M x V M H + P xVpH 

(1.5) 

P = p x V M H 
which coincides with fll.lj ) if 

H(M,P) = ±(M,n(M)) + (P,A) . (1.6) 

(Here (•, •) stands for the standard euclidean scalar product in R 3 ). The Poisson bracket ( |1.4[) has 
two Casimir functions, 

C = (M,P) and (P,P) , (1.7) 
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which are therefore integrals of motion for ( |1.1[) in involution with H{M, P) (and with any other 
function on the phase space). 

The Lagrange case of the rigid body motion (the Lagrange top, for brevity), is characterized 
by the following data: J\ = J2, which means that the body is rotationally symmetric with respect 
to the third coordinate axis), and A\ = A2 = 0, which means that the fixed point lies on the 
symmetry axis. Choosing units properly, we may assume that 

J 1 = J 2 = 1, J 3 = a, A = (0,0, 1) T . (1.8) 

The system has in this case an additional integral of motion, 

M 3 = (M,A), (1.9) 

which is also in involution with H(M,P), and assures therefore the complete integrability of the 
flow (|l.l| ). For an actual integration of this flow in terms of elliptic functions see, e.g., fGj , ||KSl , 



and for a more modern account [RM], |Au|, |CB|. 



Remarkable as it is, this result is, however, somewhat unsatisfying from the practical point of 
view. Indeed, one is usually interested in describing the motion of the top in the rest frame, which 
does not move in the physical space. It is less known that for the Lagrange top the corresponding 



equations of motion are also very nice and, actually, even somewhat simpler than (1.1): 

rh = p x a 

(1.10) 

a = m x a 

Here m = (mi,m2,ms) T € R 3 is the vector of kinetic momentum of the body, expressed in the 
rest frame, p is the unit vector along the gravity field, also expressed in the rest frame, so that it 
becomes constant: 

p = (0,0,l) T , (1.11) 

and a = (a±, 02, 03)^ € K 3 is the vector pointing from the fixed point to the center of mass, expressed 
in the rest frame. An exterior observer is mainly interested in the motion of the symmetry axis of 
the top, which is described by the vector a. 

The system ( 1.1 0| ) has several remarkable features. First of all, it does not depend explicitly 



on the anisotropy parameter a of the inertia tensor. Second, it is Hamiltonian with respect to the 
Lie-Poisson bracket of e(3): 

{mi,mj} = -e ijk m k , {mi,aj} = -e ijk a k , {ai,a,j} = 0. (1-12) 

For an arbitrary Hamilton function H(m,a), the Hamiltonian equations of motion in this bracket 
read: 

rh = V m H x m + V a H x a 

(1.13) 

d = V m H x a 
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These equations coincide with ( 1.10 ), if 



Ho(m,a) = -(to, to) + (p, a) . 



;i.i4) 



Of course, the functions 

c = (to, a) and {a, a) , (1-15) 

are Casimirs of the bracket (1.12), and therefore are integrals of motion for (1.10) in involution 
with Ho(m, a) (and with any other function on the phase space). An additional integral of motion 



in involution with Ho(m,a), assuring the complete integrability of the system (1.10), is: 



TO 3 



(m,p) . 



;i.l6) 



In the main text we give a Lagrangian derivation of equations of motion ( |1.1[ ) and fll-lCQ and an 
explanation of their Hamiltonian nature and integrability. Then we present a discrete Lagrangian 
function generating two maps approximating (LI) and ( LICD , respectively. Most beautiful is the 
discretization of ( |1.10| ): 

' m k+1 - TO fc = ep x a k , 

(1-17) 

afc+i - a k = - m k+ i x (a k + a k+1 ) . 



It is easy to see that the second equation in ( |1.17 ) can be uniquely solved for a k+ i, so that Ql.l7| ) 
defines a map (m k ,a k ) i— > (mk+i, ctfe+i), approximating, for small e, the time e shift along the 



trajectories of (1.10). This distinguishes the situation from the one in [MV] where Lagrangian 
equations led to correspondences rather than to maps. We shall demonstrate that the map (1.17) 



is Poisson with respect to the bracket ( p.. 12 ), so that the Casimirs ( 1.15| ) are integrals of motion. It 
is also obvious that ( Ljjj ) is an integral of motion. Most remarkably, this map has another integral 
of motion - an analog of the Hamiltonian: 



H £ (m,a) = ^(to,to) + {a,p) + |(o x m,p) . 



;i.!8) 



The function ( PI) is in involution with (|1.16|) , which renders the map (|1.17| ) completely integrable. 
A similar discretization for the equations of motion in the body frame ( |1.1| ) is slightly less elegant. 



2 Lagrangian mechanics on TG 
(continuous time case) 

Let li(g,g) : TG \—* K be a smooth function on the tangent bundle of the Lie group G, called 
the Lagrange function. For an arbitrary function g(t) : [to,ti] i— > G one can consider the action 
functional 

S= rL( 5 (t),5(t))dt. (2.1) 
Jt 
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A standard argument shows that the functions g(t) yielding extrema of this functional (in the class 
of variations preserving g(to) and g(t\)), satisfy with necessity the Euler-Lagrange equations. In 
local coordinates {g 1 } on G they read: 

dtWJ'dgi- [2 - 2) 

The action functional S is independent of the choice of local coordinates, and thus the Euler- 
Lagrange equations are actually coordinate independent as well. For a coordinate-free description 
in the language of differential geometry, see |Aj], [MR]. 



Introducing the quantities [] 

II = VgL € T*G , (2.3) 

one defines the Legendre transformation: 

(g,g)eTG^(g,U)eT*G. (2.4) 
If it is invertible, i.e. if g can be expressed through (g, II), then the the Legendre transformation 



of the Euler-Lagrange equations (2J2) yield a Hamiltonian system on T*G with respect to the 



standard symplectic structure on T*G and with the Hamilton function 

H(g,U) = (U,g)-L(g,g) , (2.5) 

(where, of course, g has to be expressed through (g, II)). Finally, we want to mention the Noether 
construction for deriving the existence of integrals of motion of the Euler-Lagrange equations from 
the symmetry groups of the Lagrange function. We shall formulate the simplest form of Noether's 
theorem, where Lagrangian functions are invariant under the action of one-dimensional groups. 
Let ( G j be a fixed element, and consider a one-parameter subgroup 

G«) = {e cC : c£M}cG. (2.6) 



Proposition 2.1 a) Let the Lagrange function be invariant under the action of G^ on TG 
induced by left translations on G: 

Ue«g,L ecU g)=L(g,g). (2.7) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

{R* g (VgL)X) = (R* g n,() ■ (2.8) 



For the notations from the Lie groups theory used in this and subsequent sections see Appendix A. 
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b) Let the Lagrange function be invariant under the action of on TG induced by right 
translations on G: 

Hge*,R e e U g)=L(g,g). (2.9) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(L*(V^L),C) = (^n,C). (2.10) 

Proof. Differentiate ( [2.7| ) (or Q2,9j )) with respect to c, set c = 0, and use the Euler-Lagrange 
equations. I 

For a detailed proof of the general version of the Noether theorem sec |A), [MR]. 



In practice, it turns out to be more convenient to work not with the tangent bundle TG, but 
with its trivializations G x g, which is achieved by translating the vector g dT g G into the group 
unit via left or right translations. 

2.1 Left trivialization 

Consider the trivialization map 

(g,n)eGx Q ^(g,g)eTG, (2.11) 

where 

g = L g *n o n = Lg-iJ . (2.12) 
The trivialization ( 2.11| ) of the tangent bundle TG induces the following trivialization of the cotan- 



gent bundle T*G: 

(g,M) £Gx 8 *K( 9 ,n) GT*G, (2.13) 

where 

n = L*_xM o M = L*U. (2.14) 
Denote the pull-back of the Lagrange function through 

L^(g,n)=L(g,g) , (2.15) 

so that 

L®(g,n) : GxjkI, 

We want to find differential equations satisfied by those functions (g(t),£l(t)) : [to,ti] i— > G x g 
delivering extrema of the action functional 



S«= / L^(g(t),n(t))dt 
Jt 

and such that 

O(i) = L g - Ht> g{t) . 

Admissible variations of (g(t),Q(t)) are those preserving the latter equality and the values g(to), 

g(h). 
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Proposition 2.2 The differential equations for extremals of the functional S^' read: 



where 

If the Legendre transformation 



M = ad* n ■ M + d' g L® 



9 — Lg*£l , 

M = V^L^ G g* 



(2.16) 

(2.17) 
(2.18) 



(g, n) G G x g ^ (g, M) G G x g* 

zs invertible, it turns ( [2.16 ) into a Hamiltonian form on G x g* wrat/i i/ie Hamilton function 

H(g,M) = (M,fi)-L%,n) , (2.19) 

(where, of course, f2 /ias to 6e expressed through (g,M)); the underlying invariant Poisson bracket 
on G x g* is i/ie pull-back of the standard symplectic bracket on T*G, so that for two arbitrary 
functions /^(g, M) : G x g* i— > K we have: 



{/i, h) = -«/i,Vm/ 2 ) + (d'f 2 , V M /i) + (M, [V M /i, V m / 2 ] } . 



(2.20) 



Proof. The equations of motion ( p. 16 ) may be derived by pulling back the equations ( |2.2| ) under 
the trivialization map ( 2,11| ), but it is somewhat simpler to derive them independently. To this 
end, consider the admissible variations of (g(t),Q(t)) in the form 



]{t, e) = g(t)e er >® , where V (t) : [to, h] ^ g , r/(t ) = 7?(ti) = 



and 



fi(t,e) = L ff -i (t , e) ^(t,€) = Ade^W ■ n(t) + ei^t) + 0(e 2 
= Q(t) + e^(t) + [Q(t), V (t)})+0(e 2 ). 

So, equating the variation of action to zero, we get: 

dS® 







de 



t=0 



J' 1 (( d' g L® ,n) + (V n L« , 7) + ad n • r? )) dt . 



Integrating the term with i) by parts and taking into account r/(to) = ^(^l) = 0) we come to: 

rtl «L« + ad*0-VnL«- 6 

'to 

Due to arbitariness of n(t) the following equation holds: 

d 



t ( d'L® + ad* • VnL« - ^(V n L«), 77 ) = 

Jin «t 



dt 



(V n LW) =ad*n- V n LW + d'L« . 
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It remains to notice that M denned by ( |2.14| ), ([2.3j), i.e. M = L*V^L, coincides with 
follows from ( |2,12| ). I 

We now observe, what does Noether's theorem (more exactly, its version in Proposition 2.1) 
yield under left trivialization. 

Proposition 2.3 a) Let the Lagrange function L^)(g, f2) be invariant under the action of G^ 
on G x g induced by left translations on G: 

L«(eS^)=L (0 (s^)- (2-21) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(Ad* g- 1 • V n L«, C) = (M, Adg- 1 ■ () . (2.22) 

b) Let the Lagrange function (g, Q) be invariant under the action of G^' on G x g induced 
by right translations on G: 

L«( 5 e cC ,Ade-^-ft) =L®(g,Q) . (2.23) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(V n L«C> = <M,C>. (2.24) 

We finish this subsection by discussing the reduction procedure relevant for later applications. 
Let us assume that there holds a condition somewhat stronger than ( [2.21D , namely, that the function 
is left invariant under the action of a subgroup somewhat larger than Fix an element 

C G 0, and consider the isotropy subgroup G^ of £ with respect to the adjoint action of G, i.e. 

G [c] = {h : Ad/i-C = C}cG. (2.25) 

Obviously, G^' C G^K Suppose that the Lagrange function L^(g,0) is invariant under the action 
of G^ onGxg induced by left translations on G: 

L« (hg, n) = LW (g, n) , he G [cl . (2.26) 

We want to reduce the Euler-Lagrange equations with respect to this action. As a section (G x 
we choose the set g^ x g, where g^ is the orbit of £ under the adjoint action of G: 

g c = {Adg-(, geG} Cg. (2.27) 

We define the reduced Lagrange function jC."' : g^ x g \— ► M as 

£ {l) (P,n) = L^(g,n) , where P = Ad 5 ~ 1 -C. (2.28) 

This definition is correct, because from 

p = Ad 5 r 1 -c = Ad 52 - l .c 

there follows Ad^jflf 1 • C = C, so that g 2 g± 1 € G^, and = L^(g 2 ,Q). 
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Proposition 2.4 Consider the reduction (g,Q) i— > (P, O). T/ie reduced Euler-Lagrange equations 
( pip read: 

M = ad* Jl • M + ad* P • V P £® , 

(2.29) 

P=[P,0] , 
where 

M = V n £ (,) G fl * . (2.30) 

ij i/ie Legendre transformation 

(P, fi)£g ( x B H (P, M) G 0C x s * (2.31) 

zs invertible, it turns ( p. 29 ) into a Hamiltonian system on x g*, uraia the Hamilton function 

H(P, M) = (M, Q) - £ W (P, U) , (2.32) 

where has to be expressed through (P,M); the underlying invariant Poisson structure on x g* 
is given by the following formula: 

{Pi, P 2 } = -(VpPi, [P, V M P 2 ] ) + (V P P 2 , [P, V M Pi] ) + (M, [VjifFi, V M P 2 ] ) (2.33) 

/or too arbitrary functions Pi^P, M) : g(Xg* m i. (This formula indeed defines a Poisson bracket 
on all of q x q*). 

In addition to the integral of motion ( |2.32|) , the equations of motion ( |2.29| ) always have the 
following integral of motion: 

C = (M,P) . (2.34) 



This function is a Casimir of the bracket ( |2.33 ) 



Proof is a consequence of the following formulas: 

d'gL® = ad* P • V P £« , V Q L« = Vn£ (Z) , 

which are easy to derive from the definitions, and similar formulas connecting the Lie derivatives 
of f± } 2 with the gradients of F\^- I 

2.2 Right trivialization 

All constructions in this subsection are absolutely parallel to those of the previous one, therefore 
we restrict ourselves to the formulation and omit all proofs. 
Consider the trivialization map 

SGxgM {g,g) G TG , (2.35) 

where 

g = R g „uj ^ uj = R g -i*g . (2.36) 
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This trivialization of the tangent bundle TG induces the following trivialization of the cotangent 
bundle T*G: 

(g, m) e G x q* (g, n) G T*G , (2.37) 

where 

n = Rg^m o m = R* g U . (2.38) 
The pull-back of the Lagrange function is denoted through 

lW( S)W )=%5). (2.39) 



Proposition 2.5 The differential equations for the extremals of the functional 



S^= / lS r \g(t),u{t))dt 



read: 
where 

If the Legendre transformation 



to 

rh = —ad* u • m + dgL^ , 
9 = R g *(jJ , 

m = V^LM G * . 



(2.40) 



(2.41) 



(g,Lo) £G x 2 ^ (g,m) eG x g* (2.42) 

is invertible, it turns ( [2.40 ) into a Hamiltonian form on G x g* with the Hamilton function 

%m) = (m 1 w)-LW( 9) w) ) (2.43) 

where u has to be expressed through (g,m); the underlying invariant Poisson bracket on G x g* 
is the pull-back of the standard symplectic bracket on T*G, so that for two arbitrary functions 
f 1,2(9,™) '■ G x g* 1— > R we have: 

{/l,/ 2 } = -(d g fi, V m / 2 ) + (d g f 2 ,V m fl) - (m, [V m /i, V m / 2 ] ) • (2.44) 

A version of Noether's theorem takes the following form: 

Proposition 2.6 a) Let the Lagrange function li^(g,u) be invariant under the action of G^ 
on G x g induced by left translations on G: 

(e cC 5 , Ad e cC • uj) = L< r > (g, u) . (2.45) 

TTien the following function is an integral of motion of the Euler-Lagrange equations: 

(VJLM = ^0- (2-46) 
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b) Let the Lagrange function ~L^ r \g,ui) be invariant under the action of G^ on G x g induced 
by right translations on G: 

L«( S e< )W )=LM( 9)W ). (2.47) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(Ad* g ■ vM r) ,0 = (m,Adg-(). (2.48) 

Turning to the reduction procedure, suppose that the Lagrange function l,( r \g,uj) is invariant 
under the action of on G x g induced by right translations on G: 

L W (gh, w) = LM (g, u) , he G [cl (2.49) 

We define the reduced Lagrange function : g^ x g i— > K as 

6 r) (a, u) = (g, u) , where a = Ad g ■ ( . (2.50) 



Proposition 2.7 Consider the reduction (g,uj) i— > (a,u;). T/ie reduced Euler-Lagrange equations 
(|2~40| ) read: 

m = —ad* uj ■ m — ad* a • V a £^ , 

(2.51) 

d = fa;, al , 



lodere 

7/ the Legendre transformation 



m 



V U C^ E g* . 



(2.52) 



(a,u) £ 5( x j h (a, m) £§( x g* (2.53) 

is invertible, it turns ( p. 51 ) into a Hamiltonian system on g^ x g*, wdid iae Hamilton function 

H(a,m) = (m,u) - &\a,u) , (2.54) 

where ui has to be expressed through (a,m); £/ie underlying invariant Poisson structure on g^ x g* 
is given by the following formula: 



{F 1 ,F 2 } = (V Fi, [a, V m F 2 ] } - (V a F 2 , [a, V m Fi] } - (m, [V m F 1; V m F 2 ] 



(2.55) 



/or too arbitrary functions Fi^(a, m) : g^ x g* i— > M. (This formula indeed defines a Poisson bracket 
on all of g x g*^. 

In addition to the integral of motion ( |2.54 ), the equations of motion ( p. 51 ) always have the 
following integral of motion: 

c = (m,a) . (2.56) 
This function is a Casimir of the bracket ( [2.55 ). 
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Notice that the brackets Q2.33 ) and ( 2.55 ) essentially coincide (differ only by a sign). 

Remark. For future reference notice that the elements n,wGg and M, m£j* are related via 
the formulas 

Sl = Adg~ 1 -w, (2.57) 

M = Ad* g-m. (2.58) 

3 Lagrangian mechanics on G x G 
(discrete time case) 



We now turn to the discrete time analog of the previous constructions, introduced in Q, |MV|. 
Our presentation is an adaptation of the Moser-Veselov construction for the case when the basic 
manifold is a Lie group. We shall see that almost all constructions of the previous section have 
their discrete time analogs. The only exception is the existence of the "energy" integral (|2.5| ). 

Let L(<7, g) : G x G be a smooth function, called the (discrete time) Lagrange function. For an 
arbitrary sequence {g k € G, k = ko, ko + 1, . . . , k±} one can consider the action functional 

fci-i 

S = H9k,9k+i) ■ (3.1) 
k=ko 

Obviously, the sequences {g k } delivering extrema of this functional (in the class of variations 
preserving g ko and g kl ) , satisfy with necessity the discrete Euler-Lagrange equations: 

VMg k ,g k+1 ) + V 2 Ugk-i,9k) = . (3.2) 

Here ViL(g, g) (V2L(fl r , <?)) denotes the gradient of h(g,g) with respect to the first argument g 
(resp. the second argument g). Notice that in our case, when G is a Lie group and not just a 
general smooth manifold, the equation ( p.2[ ) is written in a coordinate free form, using the intrinsic 
notions of the Lie theory. This is opposed to the continuous time case and somehow underlines the 
fundamental character of the discrete Euler-Lagrange equations. 

The equation (|3.2| ) is an implicit equation for g k +i- In general, it has more than one solution, 
and therefore defines a correspondence (multi-valued map) (g k -i,9k) l—> (9k,9k+i)- To discuss 
symplectic properties of this correspondence, one defines: 

n fc = v 2 h(g k ^ 1 , gk ) g r; k G . (3.3) 

Then ( |3.2j ) may be rewritten as the following system: 

Ilfc = -ViL(g k ,g k+ i) 

(3-4) 

Ilfc+i = V 2 L(s'fc,.S f fc+i) 

This system defines a (multivalued) map (g k , H k ) i— > (g k+ i,H k+ i) of T*G into itself. More precisely, 
the first equation in ( |3.4[) is an implicit equation for g k +i, while the second one allows for the explicit 
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and unique calculation of knowing g k and g k +i- As demonstrated in @, [MV], this map 



T*G i— > T*G is symplectic with respect to the standard symplectic structure on T*G. 

For discrete Euler-Lagrange equations there holds an analog of Noether's theorem. Again, we 
give only the simplest version thereof. 

Proposition 3.1 a) Let the Lagrange function be invariant under the action of (Jw on G x G 
induced by left translations on G: 

L(eS, e cC ff)=L( 9 ,9). (3.5) 
Then the following function is an integral of motion of the discrete Euler-Lagrange equations: 

{d 2 Ug k - 1 ,g k ),0 = (R* gk ^k,C) • (3.6) 

b) Let the Lagrange function be invariant under the action of G^' on G x G induced by right 
translations on G: 

h(ge^,ge c <)=h(g,g) . (3.7) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(-dm 9k ,g k+1 ), C) = (L* gk U k , C) ■ (3.8) 

Proof. Since both statements are proved similarly, we restrict ourselves to proof of the first one. 
To this end differentiate Q3,5|) with respect to c and set c = 0. Writing (g k ,g k+ %) for (g,g), we get: 



(diL(g k ,g k+ x),() + (^(flfc, 9k+i), C) = . 
But the discrete Euler-Lagrange equations imply that 

dih(g k ,g k+ i) = -<hli(g k -i,gk) ■ 

Hence 

(d 2 Hg k ,gk+i),C) = {d2^{g k -i,9k),Cj , 

and the statement is proved. I 

Notice that the expressions of the Noether integrals in terms of (g, IT) are exactly the same as 
in the continuous time case. 

3.1 Left trivialization 

Actually, the tangent bundle TG does not appear in the discrete time context at all. We shall 
see that the analogs of the "angular velocities" f2, to live not in T g G but in G itself. On the 
contrary, the cotangent bundle T*G still plays an important role in the discrete time theory, and 
it is still convenient to trivialize it. This subsection is devoted to the constructions related to the 
left trivialization. 
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Consider the map 

(g k ,W k ) £GxGh (g k ,g k+1 ) eGxG, (3.9) 

where 

9k+i = 9kW k Wfc = fl£ 1 g k+1 . (3.10) 

The group element W k is an analog of the left angular velocity 12 from ( |2.12| ), more precisely, it 
approximates e en . Consider also the left trivialization of the cotangent bundle T*G: 

(g k ,M k ) eGx e *H ( 9k ,U k ) G T*G , (3.11) 

where 

U k = L*^M k e> M k = L* gk U k . (3.12) 



Denote the pull-back of the Lagrange function under ( [3.9D through 

h (l) (g k ,W k )=h(g k ,g k+1 ) . (3.13) 

We want to find difference equations satisfied by the sequences {(g k ,W k ) ,k = ko,---,ki — 1} 
delivering extrema of the action functional 

and satisfying W k = g^ g k +i- Admissible variations of {(g k ,W k )} are those preserving the values 
of g ko and g kl = g kl -\W kl -i. 

Proposition 3.2 The difference equations for extremals of the functional S^ 1 ' read: 

' Ad* W k l ■ M k+1 = M k + d' g L®(g k ,W k ) , 

(3.14) 

9k+l = 9kW k , 
where 

M k = d^L<9(0 fc _i, W k -x) G q* . (3.15) 

If the "Legendre transformation" 

(g k -i, W k -i) G G x G ^ (g k , M k ) G G x g* , (3.16) 

where g k = g k -iW k -x, is invertible, then ( 3.14| ) defines a map (g k ,M k ) i— > (g k+ i, M k+ i) which is 
symplectic with respect to the Poisson bracket fl2.20p on G x q* . 

Proof. The simplest way to derive ( p. 14 ) is to pull back the equations ( |3.2j ) under the map (p.9|). 
To do this, first rewrite (|3.2| ) as 

d[t(g k ,g k+1 ) + d^L(gk-i,9k) = . (3.17) 
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We have to express these Lie derivatives in terms of (g, W). The answer is this: 

d!^{gk-ugk) = dwiPXg k -x,W k -i) , 

d[Hgk,9k+i) = d'gL®(g k ,W k ) - d w lP\g k ,W k ) . 
Indeed, let us prove, for example, the (less obvious) ( |3 . 1 9| ) . We have: 



(3.18) 
(3.19) 



(d[h(g k ,g k+1 ),rj) = — h(g k e eri , g k+l ) 



e=0 



d_ 

de 



^){g k e^,e-^W k ] 



e=0 



= (d' g L,®(g k ,W k ),ri) - {d w L®(g k ,W k ), V ) . 

It remains to substitute ( 3.18j ), ( |3.19| ) into ( 3.17Q . Taking into account that 

d w lS l \g k , Wfc) = Ad* W^ 1 • d' g \S l \g k , W k ) 

we find ( |3.14 ). Finally, notice that the notation ( |3.15 ) is consistent with the definitions (0.31), 
(|3.12| ). Indeed, from these definitions it follows: M k = d^L{gk-ii 9k), and the reference to (|1D 
finishes the proof. I 

We now observe, what does the discrete time version of the Noether theorem from Proposition 
3.1 yield under left trivialization. 

Proposition 3.3 a) Let the Lagrange function h^ l \g,W) be invariant under the action of 
on G x g induced by left translations on G: 

(e cC g, W) = L« (g, W) . (3.20) 

Then the following function is an integral of motion of the Euler-Lagrange equations: 

(Ad*^ 1 • 4dL% fc _ 1; Wfc-0,0 = (M k , Ad g^ 1 ■ . (3.21) 

b) Let the Lagrange function h^ l \g,W) be invariant under the action of G^ on G x G induced 
by right translations on G: 

lS l \ge<, e-<We<) = L«(<?, W) . (3.22) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(d^L«(ffc_ 1> W fc _i),C) = (M k ,() ■ (3.23) 

We discuss now the reduction procedure. Assume that the function is invariant under the 
action of G^ onGxG induced by left translations on G: 



lS l \hg,W) =lF>(g,W) , h€G®. 
Define the reduced Lagrange function A^ : x G i— > K as 

A^(P,W) =lS l \g,W) , where P = Adg- 1 -(. 



(3.24) 
(3.25) 
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Proposition 3.4 Consider the reduction (g,W) t— > (P, W). The reduced Euler-Lagrange equations 
fl3l4l ) read: 

Ad* W^ 1 • M fe+1 = M k + ad* P fc • V P A«(P fc , W fc ) , 

(3.26) 

P k+1 = AdW^ l -P k , 
where 

M k = d^iP^W^) G * . (3.27) 

if the "Legendre transformation" 

(P fe _!, Wjfc_i) G 8( xGk (P fc ,M fc ) G C x fl * , (3.28) 

w/iere P fc = AdW fc -\ • P fe _i, is invertible, then ( |3.26|) defme amap (P k ,M k ) i-> (Pfc + i,M& +1 ) of 
x g* which is Poisson with respect to the Poisson bracket ( |2.33 ). 
The equations of motion ( 3.26] ) always have the following integral of motion: 

C={M k ,P k ), (3.29) 

which is a Casimir function of the bracket ( 2.33| ). 

3.2 Right trivialization 

Consider the map 

(g k ,w k ) gGxGh (g k ,g k+1 ) eGxG, (3.30) 

where 

gk+i = w k g k w k = g k+ ig k 1 . (3.31) 
Consider also the right trivialization of the cotangent bundle T*G: 

Gfc, m*) £ G x g* h (g k ,U k ) G T*G , (3.32) 

where 

U k = R* ;1 m k o m fc = P* fc II fe . (3.33) 
Denote the pull-back of the Lagrange function under ( |3.30| ) through 

^ r \gk,w k ) = L(g k ,g k+1 ) . (3.34) 



Proposition 3.5 The difference equations for extremals of the functional 

fci-i 

S {r) = J2^ ir) (9k,w k ), 
feo 
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read: 

Ad* w k ■ m k+ i = m k + dglS^ (g k , w k ) , 
9k+\ = w k gk , 



(3.35) 



where 

m k = d w lS r) (g k _i,w k „i) G g* . (3.36) 

If the "Legendre transformation" 

(gk-i,Wk-l) EGxGh (g k ,m k ) € G x g* , (3.37) 

where g k = w k -x9k-x> * s invertible, then ( |3.35|) de fine a map (g k ,m k ) i— > (<?jfc+i, Tn k +i) which is 
symplectic with respect to the Poisson bracket Q2.44 ) on G x g* . 



Proof. This time the discrete Euler-Lagrange equations ( |3.2| ) are rewritten as 

diL(g k , gk+i) + d 2 L(5fe_i, g k ) = , (3.38) 

and the expressions for these Lie derivatives in terms of (g, w) read: 

dJL{gk-i,9k) = , (3.39) 

diL(g k ,g k+ i) = d g h (r \g k ,w k ) - d^V^ (g k ,w k ) = d g \S r) (g k , w k ) - Ad* w k ■ d w U r \g k ,w k ) . (3.40) 
The expression ( |3.36| ) is consistent with the definitions ( |3.3j ), ( |3.33| ), which imply that m k = 



d2^(gk-ii9k)i an d a reference to (|3.39| ) finishes the proof. I 

Proposition 3.6 a) Let the Lagrange function h^ r \g,w) be invariant under the action of G^ 
on G x G induced by left translations on G: 

L (r) (e cC g, e c( we~ c( ) = L (r) (g, w) . (3.41) 

Then the following function is an integral of motion of the Euler-Lagrange equations: 

(4l (r) (ft-i,^-i),() = (m k ,C) ■ (3.42) 

b) Let the Lagrange function h^ r \g,w) be invariant under the action of G^ on G x g induced 
by right translations on G: 

lS r \ge c ^w) =lS r) {g,w) . (3.43) 
Then the following function is an integral of motion of the Euler-Lagrange equations: 

(Ad*0 fc .d w LW(0 fe _i,w fc _ 1 ),C) = (m k , Ad g k ■ () . (3.44) 
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Finally, we turn to the reduction procedure. Assume that the function LA r ' is invariant under 
the action of G^ on G x G induced by right translations on G: 

!>) (gh, w) = I>) {g, w) , he . (3.45) 

Define the reduced Lagrange function A^ r ) : g^ x G i— > K as 

A (r) (a, w) = L (r) (g, w) , where a = Ad g ■ ( . (3.46) 



Proposition 3.7 Consider the reduction (g,w) t— > (a,w). The reduced Euler-Lagrange equations 
( 3.35| ) read: 

Ad* w k ■ m k+1 = m k - ad* a k ■ V a A^(a k ,w k ) , 



(3.47) 



a k+ i = Adw k - a k 



where 

m k = d w A^ r \a k -i,w k -i) 6 g* . (3.48) 

// i/ie "Legendre transformation" 

(ctfc-i, £9(XGh (a k ,m k ) £ §( x g* , (3.49) 

where a k = Adw k _i ■ a k _i, is invertible, then ( |3.47l ) define a map (a k ,m k ) t—> (a k+ i,m k+ i) of 



x g* which is Poisson with respect to the bracket ( 2.55| ). 



T/ie equations of motion (|3.47| ) always have the following integral of motion: 

c = (m k ,a k ) , (3.50) 



which is a Casimir of the bracket ( |2.55| ) . 

A table summarizing the unreduced and reduced Lagrangian equations of motion, both in the 
continuous and discrete time formulations, is put in Appendix B. 

4 Lagrangian formulation of the Lagrange top 

From now on we always work with the group G = SU(2), so that g = su(2), see Appendix C for 
necessary background. In particular, we identify vectors from R 3 with matrices from g, and do not 
distinguish between the vector product in R 3 and the commutator in g. We write the adjoint group 
action as a matrix conjugation, and the operators L*, R* as left and right matrix multiplication by 
<7 _1 , in accordance with ((14) and (|C.10 ). 



The following table summarizes the integrals of motion and the reductions following from the 
symmetries of Lagrange functions, in the terminology of the rigid body motion. 
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Left symmetry 
g e cp g 
(rotation about p , 
the gravity field axis) 


Right symmetry 

r A 

g i-> ge CA 
(rotation about A , 
the body symmetry axis) 


Left trivialization 
(g,U)^(g,M = g- 1 U) 
(body frame) 


(M,P), P = g~ 1 pg 


(M,A) 


Right trivialization 

0,n) ^ {g,m = iig- 1 ) 

(rest frame) 


(m,p) 


(m,a) , a = gAg' 1 



4.1 Body frame formulation 

For an arbitrary Lagrangian system on TG, whose Lagrange function may be written as 

L(g,g) = C^(P,n) , 

where Q, = g~ l g, P = g~ 1 pg, the Euler-Lagrange equations of motion take the form 

M = [M,n] + [Vp£W,P] 

P = [p,n] , 



(4.1) 



where M = VqC^ . Such systems are characterized by the condition of invariance of L(g, g) under 
the action of G^ on TG induced by left translations on G, i.e. 

L(J*g,e<*g)=L(g,g). 

The geometrical meaning of this action is the rotation around p - the symmetry axis of the gravi- 
tation field. Consider the Lagrange function of the general top: 



£«\p,n) = ^(,m,n)-(p,A), 

where J : g i— > g is a linear operator, and A £ g is a constant vector. We calculate: 

M = V n £ (I) = , VpC {l) = -A , 



(4.2) 



so that (4.1) takes the form 



where 



M = [M, n] + [P, A] , 

P = [p,n] , 
m = jn , 
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(4.3) 
(4.4) 



which is identical with ( |1.1|) . According to Proposition 2.4, this system is Hamiltonian with respect 
to the bracket ( |2,33| ), which in our case has the coordinate representation (1.4). 

The Lagrange top is distinguished by the relations fll.Sj ). They may be represented in the 
following, slightly more invariant fashion: 

M = Jtt = Q - (1 - a)(n,A)A , (4.5) 

i.e. J acts as multiplication by a in the direction of the vector A, and as the identity operator in 
the two orthogonal directions. This allows us to rewrite (|4.2[) as 



L(g, g) = (P, fi) = - (O, J2> - (O, A) 2 - (P, A) . (4.6) 



In this case the equations of motion ( |4.3| ) clearly imply that the following function is an integral of 
motion: 

C = (M, A) . 

This assures the complete integrability of the Lagrange top. 

Remark 1. It is easy to see that (|4.5| ) implies (M,A) = a(fl,A), which allows us to invert 



(4.5) immediately: 

n = M+ — - (M, A) A . (4.7) 

a 

For futher reference, we rewrite this as 

1 1 - a 

O = - M + [A, [A, Ml] . (4.8) 

a a 

This, in turn, allows us to reconstruct the motion of the frame g(t) through the motion of the 
reduced variables M(t), P(t) (actually only through M{t)). To this end one has to solve the linear 
differential equation 

g = gtt . 



Remark 2. As almost all known integrable systems, the Lagrange top has a Lax representation 



I RSTS ] , | An], the original references are [ AM |, [|R|, | RM | . It is straightforward to check the following 
Lax representation for ( [4.3| ), ( [4.5| ) with the matrices from the loop algebra sit(2)[A]: 



L(A) = [L(A),C/(A)] 



(4- 



where 



L(A) = X 2 A + AM + P , U(\) = XA + Q 



(4.10) 
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4.2 Rest frame formulation 



With the formula (| 



we can clearly rewrite the Lagrange function only in terms of u = gg 



1 1 — 

Mg,g) = £ ir \a,uj) = -{u,u) — - 



a 



{uj,a) 



(p, a) 



(4.11) 



The possibility to represent li(g,g) through u, a is equivalent to the invariance of L(g,<?) under 
the action of G^ A > on TG induced by right translations on G: 



L(ge cA ,ge cA )=L(g,g). 

The geometrical meaning of this action is the rotation around A - the symmetry axis of the top. 
The Euler-Lagrange equations of motion for such Lagrange functions read: 



m 



[Lj,m] + [a,V a £ {r) ] 



(4.12) 



to, a 



where m = V^jC^K We calculate for the Lagrange function (4.11): 

m = VujC^ = uj — (1 — a){uj, a)a , 
V a C {r) = -{l-a){u,a)u-p. 

Putting this into ( [4.12 ), we find: 

J m = [p, a] 
1 a = [m, a] 



(4.13) 



(4.14) 



which is identical with ( 1.10| ). According to Proposition 2.7, this system is Hamiltonian with respect 
to the bracket ( |2.55| ), whose coordinate representation coincides with ( |1.12[ ). 

Remark 1. It follows from ( |4 . 1 3| ) that (m,a) = a(uj,a), so that ( |4.13j ) can be easily inverted: 



1 — a , 

uj = m H {m, a) a . 

a 



(4.15) 



Recall that c = (m,a) is a Casimir function of the underlying invariant Poisson bracket ( 1.12] ). Now 
the latter formula allows us to reconstruct the frame evolution from the evolution of the reduced 
variables (m, a) via integration of the linear differential equation 



9 = ug 



Remark 2. It turns out to be possible to derive from ( 4.14| ) a closed second order differential 
equation for a. Indeed, take the vector product of the second equation in Ql.lCp in order to obtain 



m 



a x a + ca 



(4.16) 
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Substituting this into the first equation in ( p.. 10 ), we find: 

axa + ca = pxa 



(4.17) 



Remark 3. The Lax representations for (4.14) is, of course, gauge equivalent to the one for 
the body frame formulation, but is slightly simpler than the latter Q, [RSTS]. It reads: 



i(x) = [e(x), u (\)}, 



(4.18) 



with the matrices 



£(X) = X 2 a + Xm + p , u(X) = Xa . 



(4.19) 



In Sect. 6 we indicate how this Lax representation can be derived from the zero curvature represen- 
tation of the so called Heisenberg magnetic. 



5 Discrete time Lagrange top 

We now give (in an ad hoc manner) the discrete Lagrange function which is claimed to lead to a 
suitable discretization of the Lagrange top. The motivation for the choice of this function comes 
from the geometry of curves and will be given in the next section. Unlike the continuous time case, 
we start with the rest frame formulation. 

5.1 Rest frame formulation 

Consider 

H9k,9k+i) = ^ {r) (ak,w k ) = -^logtr(wfc) - — - °^ log (l + (a k , w k a k w^ 1 )^ - e(p,a k ) , (5.1) 

where a k , w k are defined as in Sect. 3.2: w k = g k +\g k , a k = g k Ag~^ x . Notice that {ak^WkakW^ 1 ) 
in (|5.1|) is nothing but (a k , a k+ \). To see that the function (|5.1|) indeed gives a proper discretization 



of (4.11), we shall need the following simple lemma. 

Lemma 5.1 Let w(e) = 1 + ecu + 0(e 2 ) G SU(2) be a smooth curve, uj G su{2). Then 



tr(w(e)) = 2-j(cu,u;) + 0(e 3 ) . (5.2) 
For an arbitary a G su{2): 

{a,w(e)aw' 1 (e)) = {a, a) - —((a,a){u),u) - (a,w) 2 J + 0(e 3 ) . (5.3) 
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Proof. Let w = 1 + ecu + e 2 v + 0(e 3 ). Then from ww* = 1 we get: 

1 



Hence 



v + v* + ujco* = => v = ^uj 2 + v 1 , v\ G su(2) . (5.4) 



tr( v ) = ^tr(u; 2 ) = -- (w.w), 



waw* = a + e[o;, a] + — [a;, [a;, a]] + e 2 [vi, a] + 0(e 3 ) 



which proves (5.2). Similarly, we derive from (^ 

2~ 

which implies Q5.3j). I 

With the help of this lemma we immediately see that, if w = 1 + ecu + 0(e 2 ), then, up to an 
additive constant, 

A {r) (a,w) = eC {r) (a,Lo) + 0(e 2 ) , 
where £( r \a,oj) is the Lagrange function ( |4.11| ) of the Lagrange top. 



Theorem 5.2 The Euler-Lagrange equations of motion for the Lagrange function (5T) are equiv- 
alent to the following system: 

m k+ i = m k + e[p,a k ] , 

(5.5) 

flfc+i = a k + - [m k+ i,a k + a k+ i] . 
The second equation of motion can be uniquely solved for a k +\ : 

a k+ i = (1 + em k+1 )a k (l + £m w ) _1 . (5.6) 

The map (m k , a k ) \— > (mjt-fi, afc+i) is Poisson with respect to the bracket ( |1.12j ) and /ias two integrals 
in involution assuring its complete integrability: (m,p) and 

H E (m,a) = -(m,m) + (a,p) + |([a,m],p) . (5.7) 

Proof. According to Proposition 3.7, the Euler-Lagrange equations of motion have the form: 

w^m^Wk = m k + [a k , V a A^ (a k , w k )\ , 

(5.8) 

a k+1 = w k a k w k 1 , 
where 

m k+ i = d w A^(a k ,w k ) . (5.9) 
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To calculate the derivatives of A^ r \ we use the following formulas: 



d w tr(w k ) = --$s(w k ) , d w {a k ,w k a k w k 1 ) = [a k+1 ,a k ] . 

^a{a k ,w k a k w^ 1 ) = a k+1 + w k ~ 1 a k w k . 
Indeed, the first one of these expressions follows from: 



(5.10) 
(5.11) 



(d w ti(w k ),rj) = ^-tr(e ev w k ) 



e=0 



tr(r]w k ) = tr(r/3f(u; fc )) 



1 



6=0 



To prove the second one, proceed similarly: 

{d w {a k ,w k a k w k ~ 1 ),r]) = ^{a k , e tr] w k a k wl x eT^) 
Finally, as for the third expression, we have: 

(V a (a k ,w k a k w k ~ 1 ),r]) = — {a k + erj, w k (a k + erfw^ 1 ) 



{a k ,[i],a k+1 }) = ([a k+1 ,a k ],rj) . 



e=0 



(a k+1 +w k 1 a k w k ,7]) . 



With the help of fl5lc| ), ( ^TT|) we find the following expressions: 

2a %(w k ) 2(1 -a) [a k+1 ,a k ] 
m k+ i = — 



£ tv(w k ) 



e l + (a k ,a k 



+1/ 



(5.12) 



and 



w, 



~ l m k+ iw k - [a k , V a A {r \a k ,w k )] = — 



2a Q(w k ) 2(1 -a) [a k+1 ,a k ] 



£ tr(w k ) £ l + (a k ,a k+1 ) 

m k+ i + £[a k ,p] . 



+ s[ a k,P] 



(5.13) 



Comparing the latter formula with the first equation of motion in (5.8), we find that it can be 
rewritten as 

m k+ i +£[a k ,p] = m k , 

which is equivalent to the first equation of motion in (|5.5|). 

To derive the second one, rewrite the second equation in (|5.8D as 



= a k+1 w k - w k a k = ^t(w k )(a k+ i - a k ) + a k+ i^s(w k ) - Q(w k )a k 
= ^{w k ){a k+ i - a k ) + -[a k+1 + a k ,%(w k )] 

(we used Lemma C.3 and the equality (a k+ \ ,Q(w k )) = {a k ,^s(w k )) which follows from the same 
equation a k+ \w k = w k a k we started with). So, the second equation in (|5,8|) is equivalent to 



Ofc+i - a k 



tr(w k ) 



Ofc+i + a k 



(5.14) 
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On the other hand, for any two unit vectors a k , a k +i with a k +\ + at ^ we have: 

[«fc+l, ojfe] 



1 + (a fc ,a fc 



+1/ 



(5.15) 



Comparing ( |5.14 ), ( 5.15| ) with ( |5.12 ), we find the second equation of motion in (|5.5|). 

Next, we want to show how the second equation of motion in ( |5.5D can be solved for a k +i- This 
equation implies (afe+i, rrik+i) = (a k , m k+ i), so that, according to Lemma C.3, it can be rewritten 
as 

+ zak+imk+i = Ofc + £ m k+i a k , 

which is clearly equivalent to d5 



The Poisson properties of the map (5J3) are assured by Proposition 3.7. 

It remains to demonstrate that the function (|5.7| ) is indeed an integral of motion. This is done 
by the following derivation: 



H £ {m k+ i, a,k+i) 



-{m k+1 ,m k +i) + (a k+ i + | [a k+ i,m k+1 ] ,p) 
-{m k+1 ,m k+ x) + (a k - ~ [a k ,m k+1 ] ,p) 
-{m k+1 ,m k+ x - e[p,a k ] ) + (a k ,p) 
-{m k + e\p,a k ),m k ) + (a k ,p) = H £ (m k ,a k ) . 



The theorem is proved. I 

Remark. The equations of motion (|5.5|), being written entirely in terms of elements of the Lie 



algebra su(2), are clearly equivalent to the equations of motion ( 1.17 ), which are written in terms 
of vectors from M 3 . The situation with ( |5.6| ) is slightly different. Indeed, it corresponds to the 
following formula in M 3 : 

1 + em k+1 /2 

a k+l — Vfc+lOfc — Z 777 a k , 

1 - em k+1 /2 

where the orthogonal matrix Qfc+i € SO (3) is constructed out of the skew-symmetric matrix 
m k+ i € so (3) which corresponds to the vector m k +i 6 R 3 according to the following rule: 



m 



(mi,m 2 ,m 3 ) T G 



m 



—ms 
1TL3 
— m,2 mi 




e so(3) . 



Just as in the continuous time case, it is possible to derive a closed second order difference 
equation for the motion of the body axis a k . 
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Proposition 5.3 The sequence of a k satisfies the following equation: 

afcx — r + t— 7 r + ec 7—7 r -7—7 \ = e P xa k, (5-16) 

VI + {a k ,a k+ i} l + (a fc _i,a fc )/ Vl + (afc,afc+i) 1 + (a fe _i,a fc )/ 

where c = {m k ,a k ) is an integral of motion. 

Proof. Take a vector product of the second equation of motion in ( [1.17 ) by a k+ \ + a k . Taking 
into account that {m k+ i, a/c+i) = {m k+ i,a k ) = c, we find: 

2a k x afc+i = emfc + i(l + (a k , a k+ i)) - ec(a k+1 + a k ) , 

or 

2 a k x a fc+1 afc + i + 

w fc+ i = -— \+ c T~r7 \- ( 5 - 17 ) 

el + {a k ,a k+1 ) l + {a k ,a k+1 ) 

Plugging this into the first equation of motion in ( 1.17| ), we arrive at ( |5.16j ). I 

Further, we demonstrate how to reconstruct the "angular velocity" w k (and therefore the motion 
of the frame g k ) from the evolution of the reduced variables (a k ,m k ). 

Proposition 5.4 The discrete time evolution of the frame g k can be determined from the linear 
difference equation 

9k+\ = w k g k , (5.18) 

where w k are given by 

^ = t£(p) (1 + £6)) (519) 

where 

I- a a k+ i + a k 2 a k x a k+ i c a k+1 + a k 

Kk = m k+ i + c — — - = -—— + __^_^ (5.2O) 

a l + {a k ,a k+1 ) el + (a k ,a k+1 ) a 1 + (a k , a k+1 ) 

and 



, , n 2 L 1 + {a k ,a k+1 ) 

tT{Wk) = I g2 = V 2 l + e*<*/4a* ■ (5 ' 21) 

y 1 + 

Proof. We combine ( |5.12 ) with ( |5.17| ) in order to derive the formula 

2 %(wk) = £ , 
tr(w k ) 

with the expressions for £ k given in ( 5.20] ). Now the reference to Lemma C.2 finishes the proof. I 
Finally, we give a Lax representation for the map (|5.5|). 
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Theorem 5.5 The map (|5.5| ) has the following Lax representation: 

4+1 (A) = u^(X)£ k (X)u k (X) , (5.22) 

w;zf/j the matrices 

4(A) = A 2 ^a Jfe + |[o Jfc ,m fc ] + jpj + Am fc +p, u fc (A) = 1 + eAa fc . (5.23) 

Proof - a direct verification. I 

In the next section we present a derivation of this Lax representation from the one for the so 
called lattice Heisenberg magnetic. 

5.2 Moving frame formulation 

Note that the discrete Lagrange function ( |Q| ) may be also expressed in terms of P k = g k 1 pg k , 

w k = 9k l gk+i- 

Ug k ,g k+ i) = A {l] (P k ,W k ) = --^logtr(iy fc )- 2(1 ~ Q) log (l+(A,W^AW k ))-e(P kl A) . (5.24) 

Since W k = 1 + sQ + 0(e 2 ), we can apply Lemma 5.1 to see that 

A {l) (Pk,W k )=e£^(P,n) + 0(e 2 ), 

where C^ l \P, f2) is the Lagrange function ( [4.6[) of the continuous time Lagrange top. Now, one can 
derive all results concerning the discrete time Lagrange top in the body frame from the ones in the 
rest frame by performing the change of frames so that 

M k = g k X m k g k , P k = g]~ pg k , A = g^ l a k g k . 



Theorem 5.6 The Euler-Lagrange equations for the Lagrange function ( |5.24| ) are equivalent to 
the following system: 

M k+l = W^ 1 (Af* + e[P k ,A] ) W k , 



(5.25) 



Pk+i = W k l P k W k 



where the "angular velocity" W k is determined by the "angular momentum" M k+ \ via the following 
formula and Lemma C.2: 



%(W k ) e 2(1 -a) A,(l + eM k+1 )- 1 A{l+eM k+1 

7 /tt 7~ N = — M k + l H 7 

tr (Wfc) « <* 1 + (A, (1 + eM k+l )- l A{l + eM k+l 

= £ (-M k+1 + ^^[A,[A,M k+1 ]])+0(e 2 ). 



(5.26) 
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The map ( 5.25[) , ( |5.26 ) is Poisson with respect to the Poisson bracket ( |1.4| ) and has two integrals 
in involution assuring its complete integr ability: (M, A) and 



H e (M,P) = ~<M,M) + (P, A) + £ -{[M,P\,A) . 



(5.27) 



Remark. It might be preferable to express W k through (M k ,P k ) rather than through M k+ \ 



(in particular, this is necessary in order to demonstrate that the map (M k ,P k ) {M k+ \ , Pk+i ) is 
well defined). The corresponding expression reads: 



o 3(W fc ) 
tr(W fc ) 

W k AWr ] 



-(M k + e[P k ,A]) 
a 



2(1 



a 



[A, W k AW^] 



a 



i\ ' 



1 + (A, W k AW k , 
l + eM k + e 2 {P k ,A]) A(l + eM k + e 2 [P k , A])' 1 . 



(5.28) 



(5.29) 

We see that the resulting formula is similar to ( 5.26| ), but its right-hand side depends not only 
on M k but also on P k (though this latter dependence appears only in 0(e 2 ) terms). According to 
Lemma C.2, both versions allow for the reconstruction of the evolution of the frame g k from the 
evolution of the reduced variables (M k ,P k ), anyway. 

We close this section with a Lax representation for the map (5.25), ( 5.26| ). 
Theorem 5.7 The map ( 5.25| ), ( |5.26| ) has the following Lax representation: 

L k+1 (X) = C/ fe - 1 (A)L A; (A)[/ fc (A) , 



with the matrices 



L k (X)=X 2 {A+ e -[A,M k ] + £ -P k \+XM k + P k , 



[4(A) = (1 + eXA)W k 



(5.30) 



(5.31) 



Proof - a direct verification. 



6 Motivation: Lagrange top and elastic curves 



The Lagrange function (5.1) was found using an analogy between the Lagrange top and the elastic 
curves as a heuristic tool. The present section is devoted to an exposition of the corresponding 
interrelations. 

Let 7 : [0,/] i— > M 3 be a smooth curve parametrized by the arclength x G [0,1]. Defining 
the tangent vector T : [0,1] i— > R 3 as T(x) = r y'(x), the characteristic property of the arclength 
parametrization may be expressed as 

|T(x)| = l, (6.1) 
where | • | stands for the euclidean norm. The curvature of the curve 7 is defined as 

k(x) = \T'(x)\. (6.2) 
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Definition 6.1 Q, |LSf| . A classical elastic curve (Bernoulli's elastica) is a curve delivering an 
extremum to the functional 

I k 2 (x)dx , (6.3) 
Jo 

the admissible variations of the curve are those preserving 7(0) and more precisely, those 

preserving 7(Z) — 7(0) = f T(x)dx. 

Introducing the Lagrange multipliers p G M 3 corresponding to this constraint, we come to the 
functional 

j (\T'(x)\ 2 - 2{p,T(x)f)dx . (6.4) 

Identifying the arclength parameter x with the time t, this functional becomes (twice) the action 
functional for the spherical pendulum. So, classical elasticae are in a one-to-one correspondence 
with the motions of the spherical pendulum. 

A generalization of these notions to elastic rods (which physically means that they can be 
twisted) requires the curves to be framed, i.e. to carry an orthonormal frame <5(x) = (T(x), N(x), 
B{x)) in each point. In other words, a framed curve is a map <3? : [0,1] >— > {frames}. The curve 
itself is then defined by integration: 7(3;) = T(y)dy. The following quantities are attributes of a 
framed curve: the geodesic curvature 

h(x) = (T'(x),N(x)) , (6.5) 

the normal curvature 

k 2 (x) = (T'(x),B(x)) , (6.6) 

and the torsion 

r(z) = (N'{x),B{x)) . (6.7) 
Obviously, one has: k 2 (x) = kf(x) + /c|(x). 

Definition 6.2 D, fLSfl . An elastic rod (Kirchhoff's elastica) is a framed curve delivering an 
extremum to the functional 

(k 2 {x) + aT 2 {x)}dx (6.8) 

with some a / 0. The admissible variations of the curve preserve ^(O), <£(/), and j(l) — 7(0) = 
Jo T{x)dx. 

The first term in ( ]6.S{ ) corresponds to the bending energy, the second one corresponds to the twist 
energy. 

We shall identify M 3 with su(2), as described in Appendix C, and the frames with elements of 
€ SU(2), according to the following prescription: 

T = $~ 1 e 3 $, N = ^ 1 e 1 <^, B = $ _1 e 2 $ . (6.9) 
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Then, denoting 

n = , u = , (6.10) 

we find: 

k x = (u,B) = (0,e 2 ) =n 2 , fe 3 = -(a; 1 JN0 = -(n s ei) =-fli, (6.11) 

T = (u,T) = (n,e 3 ) = n 3 . (6.12) 

So, the variational problem for elastic rods may be formulated as follows: find <J>(x) : [0, 1] h- > 5C7 (2) 
delivering an extremum of the functional 

y' + nl(a;) + an§(s) - 2 (p, T(a?)>) cte , (6.13) 

where p is an (^-independent) Lagrange multiplier coming from the condition of fixed j(l) — 
7(0) = JqT(x)cIx. Identifying the arclength parameter x with the time t, <5(x) = <7 _1 (i), so that 
Q(x) = -^'(x)^- 1 ^) = g-^git) = J2(t), and T(x) = $~ 1 (x)e 3 $(x) = 5(t)e 3 ^ 1 (t) = o(t), we 
see that the functional ( 6.13j ) coincides with (twice) the action functional for the Lagrange top. 
This proves the 

Proposition 6.3 (Kirchhoff's kinetic analogy, J]]). The frames of arclength parametrized elastic 
rods are in a one-to-one correspondence with the motions of the Lagrange top. 

Actually, we use another characterization of the elastic rods. From the Euler-Lagrange equa- 
tions it follows: 



Proposition 6.4 The torsion r along the extremals of the functional ( 6.13 ) is constant, and the 
tangent vector T{x) satisfies the following second-order differential equation: 

T xT" + cT' = pxT , (6.14) 

where c = ar. Conversely, each solution T(x) of ( |6.14j ) corresponds to a curve j(x), which, being 
equipped with a frame with constant torsion t, delivers an extremum to the functional ( |6,13j ) with 
a = c/t. 

Equation ( gig ) is fl4T7|) in new notations. The latter differential equation allows the following 
interpretation. Consider the so called Heisenberg flow. It is defined by the differential equation 

T t = T xT" , (6.15) 

and describes the evolution of a curve in the binormal direction with the velocity equal to the 
curvature. Here the "time" t has nothing in common with the time t of the Lagrange top, which 
is, remember, identified with x. It is easy to see that the flow on curves defined by the vector field 
T x = T' (a reparametrization of a curve) commutes with the Heisenberg flow ( |6.15| ). Using this 
fact, we can integrate fl6.15|) once in order to find 

7t = 7 'x 7 " = TxT'. (6.16) 

(The reparametrization flow, once integrated, takes the form j x = 7' = T). Now we can formulate 
the following fundamental statement. 
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Theorem 6.5 [ Haf] , [jCq] . Let $ : [0, Z] i— > SU{2) be the frame of an elastic rod, and 7 : [0, 1] 1— > 
sn(2) i/ie corresponding curve with the tangent vector T = 7' : [0,/] 1— > sn(2). T/ien £/ie evolution 
of 7 under the Heisenberg flow ( |6,16| ) is a rigid screw-motion, and the evolution of T under the 
Heisenberg flow ( |6,15| ) is a rigid rotation. Conversely, if the evolution of T is a rigid rotation, then 
T can be lifted to a frame <i? of an elastic rod. 

The first statement of the theorem follows from ( |6.14 ). The left-hand side of ( 6.14| ) can be 
interpreted as the vector field on curves, corresponding to a linear combination of the Heisenberg 
flow and the reparametrization: 

T t + cT x =pxT . 

Integrated once, this equation yields a rigid screw motion for the curve 7: 

7t + C7x = p x 7 + q , 

where q G su(2) is a fixed vector. The converse statement follows from Proposition 6.4. 

By the way, this theorem allows to find a Lax representation for the equation ( |6.14 ), and 
therefore for the Lagrange top, starting from the well-known Lax representation for the Heisenberg 
flow. 



Proposition 6.6 The equation (6.14) is equivalent to the Lax equation 

e x (\) = [e(x),u(x)} (6.17) 

with the matrices 

e(X) = X 2 T + X{TxT x + cT)+p, u(X) = XT. (6.18) 

Proof. Indeed, the Heisenberg flow ( |6.15| ) is equivalent to the following matrix equation ("zero 
curvature representation", [|FT| ): 

ut — v x + [v, u] = , 
where u,v £ su(2)[X] are the following matrices: 

u = XT , v = X 2 T + XT xT x . 



Now it is easy to derive that the equation ( 6.14 ), rewritten as 

T t + cT x = [p, T] , 

is equivalent to ( 6.17| ) with I = v + cu + p. I 

Remembering that in the Kirchhoff 's kinetic analogy x is identified with t, T is identified with 
a, and recalling the formula (|4.16| ), we recover the Lax representation of the Lagrange top in the 
rest frame given in ( 4.18[ ), ( 4.19| ). 

Theorem 6.5 is also a departure point for discretizing elastic curves and, therefore, the Lagrange 
top 0. A discrete arc-length parametrized curve is a sequence 7 : Z 1— > R 3 with the property 
|Tfc| = 1, where Tk = jk ~ 7fc-l- Correspondingly, discrete framed curves are the sequences of 
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orthonormal frames & k , such that T k = <& k ~ 1 e3<& k 0. As before, we identify R 3 with su(2), and 
the space of orthonormal frames with SU(2). The curve 7 can be reconstructed by applying the 
summation operation to the sequence T. 

A discretization of the Heisenberg flow is well known [3kl|, [FT], see also [DS[ for geometric 
interpretation of discrete flow. It reads: 



(Tk)t 



2 T k x T k+ i 



2 7 1 fc _ 1 x T k 



l + (T k ,T k+1 ) l + (T k _ u T k ) 
A commuting flow approximating T x = T' is given by: 

T k + T k+ i Tfc_i + T k 



{Tk)x ' l + (T k ,T k+1 ) 
Once "integrated", this gives the flows on 7^: 

T k + T k+ i 



l + (T k _ u T k ) 



{lk)x 



1 + (Tk, T k+ \) 



(lk)t 



2T k x T k+ i 
1 + (T k ,T k+ i) 



(6.19) 



(6.20) 



(6.21) 



Now we accept the following discrete version of Theorem 6.5 as a definition of discrete elastic 
rods. 

Definition 6.7 A discrete elastic rod is a framed curve for which the evolution of '7^. under a linear 
combination of flows (7&)t + c(^ k ) x with some c is a rigid skrew-motion, so that the evolution of 
T k under the flow (T k ) t + c(T k ) x is a rigid rotation. 

In other words, the sequence T k satisfies the following second order difference equation: 



T k x 



+ 



21V 



l + (T k ,T k+1 ) l + (T k _ x ,T k ) 



+ c 



T k+ \ + T k 



T k + Tfc_i 



l + (T k ,T k+1 ) l + (T k _ u T k 



P xT k . (6.22) 



(This is the equation ( 5. 16| ) with e = 1 in new notations). We can immediately find the Lax 
representation for the difference equation (6.22). 



Proposition 6.8 The equation ( 6.22 ) is equivalent to the Lax equation 

£ k+1 (X)=u k 1 (X)£ k (X)u k (X) 

with the matrices 

(1 + A 2 /4) 4(A) = (A 2 + cA) Tk ,t T t~ 1 ; + (2A - cA 2 /2) ^ X Tfe 



l + (r fcj T fc _i) l + (T fc ,T fc _i) 

« fc (A) = 1 + XT k . 



(6.23) 



+ (l + A 2 /4)p, (6.24) 



(6.25) 



2 Note that the frames as well as the tangent vectors Tk, are attached to the edges [7fc-i,7fc] of the discrete 
curve 7 
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Proof. It is well known (see fFlj| ) that the flows ( 6.19| ), ( |6.20| ) allow the following "discrete zero 
curvature representations" : 

( \ (i) (i) ( \ (o) (o) 

respectively, with the matrices u k as in ( 6.25| ) and 

/I i \2/,i\ (1) \2 Tk+Tk-l . T k -i x T k 

(1 + A 2 /4)4 0) =A Tk+T ^ r » xTt 



l + (T fc ,T fc _i) 2 l + (T fc ,T fe _i) ' 
Now it is easy to see that the equation ( |6.22| ), rewritten as 

(r fc )t + c(T k ) x = \p,T k ] , 

is equivalent to u k £ k+ i = £ k u k with 

e (!) , (°) , 
4 = v k ' + cv k + p , 

which coincides with (6.24). I 

To establish a link with the discrete time Lagrange top, recall that the formula ( |5.17 ) in our 
new notations reads: 

Tfc-i x T k T k + T k _i 
mrik = 2 ; r + c ; r , 

which implies also 

^ , 1 „ Tfc + T fe _i c T fc _i x T fc 

Tfc + -x T k x m fc 



2 i + <r fcj r fc _i) 2 i + (r fcj r fc _i) ' 

Hence we can write: 

(l + A 2 /4)4 = A 2 (T fc + ir fe xm fc + ip) + Xm k + p , 

which coincides with (|5.23|) up to a nonessential constant factor. 

It remains to find a variational problem generating the equations of motion ( |6.22| ). But the 
calculations of Sect. 5 show that this task is solved by the functional (|5.1j). This gives the following 
alternative definition of discrete elastic rods. 

Definition 6.9 A discrete elastic rod is a discrete framed curve given by a finite sequence 
G su(2) delivering an extremum to the functional 

N-l N 

£ (-4alogtr(cI>- 1 <D fc+1 )-2(l-a)log(l + (r fc ,r fc+1 ))) -]T(p,T fc > (6.26) 

k=l k=l 

with some a ^ 0. The admissible variations of the curve preserve $i ; $at ; and 7jv — 70 = 
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The equivalence of Definitions 6.7 and 6.9 is the basic new result of this section. It is a geometric 
counterpart and a motivation for the considerations of Sect. 5. 

We want to close this section by giving discretizations of geometrical notions like curvature and 
torsion. Notice that the functional ( |6.26 ) naturally splits into two parts, one independent on a and 
one proportional to a. Accordingly, we declare 

- 2 J2 lo g(! + ( T k,T k+1 )) = 2 ^ log (l + hi) + const (6.27) 
k k ^ ' 

as a discretization of the "bending energy" \ Jq k 2 (x)dx, and 

Y, ( - 4 log tr^+x^- 1 ) + 21og(l + (Tfc.Tfc+i))) = 2 ]Tlog (l + ^r fc 2 ] + const (6.28) 



as a discretization of the "twist energy" 1 /QT 2 (:E)dx. Here we define the "discrete curvature" 
at the vertex 7^ by 

1 2 

1 + 7^ = ; ~ .rp nr T <= fcfc = 2tan(^ fc /2), 

where ipk is the angle between the vectors T k and T k +i- Notice that the k k depends not on the 
whole frame, but on the tangent vectots T k only, so that it makes sense also for non-framed curves. 
The "discrete torsion" r k at the vertex 7^ is defined by 

1 2 _ 2(i + (r fc ,r fc+1 )) o / s($ fc+1 ^ 1 ) \ 

+ 4 Tfc " (tr(d> fc+1 ^) )2 Tfc " 2 \tr($ fe+1 ^)' e 7 • 

The last formula will be commented on immediately. Let us demonstrate that, in a complete analogy 
with the continuous case, the discrete torsion is constant along the extremals of the functional fl6.26|) . 
Denoting for a moment 

= ( -1 S ) e sc, < 2 > • 



we find: 
and also 

so that 



1 + {T k ,T k+1 ) = 1 - 2tr($ fc+1 ^ 1 e 3 $ fc ^ 1 e 3 ) = 2|a| 2 , 

»(a) = itr^fc+i^ 1 ) , 9(a) = tr^+^ea) = -^m+i^ 1 ), e 3 ) , 

r k = 2 §4 = -2 / , -\ = -2 / ^^±1) , T fc+1 \ . (6.29) 

9(a) \tr(d> fc+1 $- 1 )' 7 \tr($ fe - 1 ^ +1 )' * + 7 1 ' 

Comparing this with fl5.12|) (remember, we set e = 1 and identified a k with T k and w k with <&~^\<& k ), 
we see that 

T fc = c/a , 
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Figure 1: Evolution of the axis of the discrete spinning top 



where c = {rn^i, 2fc.fi) is an integral of motion of the Euler-Lagrange equations (a Casimir function 
of the e(3) Lie-Poisson bracket). This corresponds literally to the continuous case. 

Remark. The case a = corresponds to discrete elastic curves 7 : Z 1— > R 3 . The tangent 
vectors T : Z \— > S 2 , Tfc = 7^. — 7fc_i, define a trajectory of the discrete time spherical pendulum. 
Its Lagrange function is obtained, as in the continuous time case, from the bending energy (|6.27 ), 
upon introducing the Lagrange multiplier p. Notice that the Lagrange function of the discrete time 
spherical pendulum is defined on S 2 x S 2 . 



7 Visualisation 

After the theory has been developed, it is tempting to look at the spinning of the discrete time 
Lagrange top. Fortunately, in the computer era, a discrete time top is even simpler to simulate than 
a classical one. Indeed, as it is shown in Theorem 5.2, the Poisson map (?7tfc,afc) 1— > (mk+i,o,k+l) 
is well defined and can be easily iterated. The vectors having been computed, Proposition 5.4 
provides us with the evolution of the frame g^, which describes the rotation of the top completely. 
So, given (mo,Oo), the rotation of the top is determined uniquely. Due to ( |5.17| ) one can take 
two consecutive positions (ao,Oi) of the axis as the initial conditions as well. Fig.l demonstrates 
a typical discrete time precession of the axis. Compare this with the classical continuous time 



pictures in |]KS| , |A|. 

The motion of the discrete time Lagrange top can be viewed using a web-browser. The Java- 
applet has been written by Ulrich Heller and can be found on the web page 

\protect\vrule widthOpt\protect\href {http : //www-sf b288 . math . tu-berlin . de\string~bobenko/bobenk< 

The applet presents an animated spinning top described by the formulas of the present paper. 
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8 Conclusion 



We took an opportunity of elaborating an integrable discretization of the Lagrange top to study 
in a considerable detail the general theory of discrete time Lagrangian mechanics on Lie groups. 
We consider this theory as an important source of symplectic and, more general, Poisson maps. 
Moreover, from some points of view the variational (Lagrangian) structure is even more fundamental 
and important than the Poisson (Hamiltonian) one (cf. ]HMR| ], |MPS |, where a similar viewpoint 
is represented). 

It is somewhat astonishing that this construction is able to produce integrable discrete time 
systems, since integrability is not built in it a priori. Nevertheless, we extend the Moser-Veselov's 
list fVf , [MV] of integrable discrete time Lagrangian systems with a new item, namely, an integrable 
discrete time Lagrange top. It seems that this list may be further continued. 

In finding this new discrete time mechanical system an analogy with some differential-geometric 
notions was very instructive. Also these interrelations between integrable differential geometry and 
integrable mechanics, both continuous and discrete, deserve to be studied further. 

Let us mention also some more concrete problems connected with this work. First of all, the 
discrete time Lax representations found here call for being understood both from the r-matrix 
point of view [RSTS], (S| and from the point of view of matrix factorizations [ MV| (unfortunately, 
these two schemes, being in principle closely related, still could not be merged into a unified one). 
Futher, the discrete time dynamics should be integrated in terms of elliptic functions. The methods 
of the finite-gap theory will be useful here [ RM |. Finally, it would be important to elaborate a 
variational interpretation of different integrable discretizations of the Euler top found in [BLS]. 



A Notations 

We fix here some notations and definitions used throughout the paper. 

Let G be a Lie group with the Lie algebra g, and let q* be a dual vector space to g. We identify 
and g* with the tangent space and the cotangent space to G in the group unity, respectively: 

= T e G , 0* = TIG . 

The pairing between the cotangent and the tangent spaces T*G and T g G in an arbitrary point g G G 
is denoted by (•,•). The left and right translations in the group are the maps L g ,R g : G i— > G 
defined by 

L g h = gh , R g h = hg V7i € G , 
and L gif , R gjf stand for the differentials of these maps: 

L g * : ThG i-> T g hG , R g * : T^G i-> Th g G . 

We denote by 

Adg = Lg^Rg-i* : ^ 
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the adjoint action of the Lie group G on its Lie algebra g = T e G. The linear operators 

L* g : T* gh G » T* h G , R* g : T* hg G » T* h G 
are conjugated to L gif , R gif , respectively, via the pairing (■, •): 

(L* g Z,n) = (Z,L g *n) for UT^G, rj € T h G , 

{R* g Z,ri) = {Z,Rg.ri) ^ £ € Tf ig G , ^T^G. 
The coadjoint action of the group 

Ad 5 = L 9 i? g -i : i-> 
is conjugated to Adg via the pairing (-, •): 

(Ad^-^r?) = (£,Ad 5 -r/) for £ £ 0* , r] E q . 
The differentials of Ad g and of Ad* g with respect to g in the group unity e are the operators 

ad r\ : i— > and ad* 77 : 0* 1— > 0* , 
respectively, also conjugated via the pairing (•, •): 

<adV£,0 = (£,ad77-C> V£ € 0* , (Gg. 
The action of ad is given by applying the Lie bracket in 0: 

ad»/-C=fo,C], VCG0. 

Finally, we shall need the notion of gradients of functions on vector spaces and on manifolds. 
If X is a vector space, and / : X 1— > R is a smooth function, then the gradient V/ : A? 1— > X* is 
defined via the formula 



(Vf(x),y)= j e f(x + ey) 



VyeX . 

e=0 



Similarly, for a function / : G <—> R on a smooth manifold G its gradient V/ : G i-> T*G is defined 
in the following way: for an arbitrary 5 G T 5 G let g(e) be a curve in G through g(0) = g with the 
tangent vector g(0) = g. Then 



(Vf(g),g)= jjm) 



e=0 



If G is a Lie group, then two convenient ways to define a curve in G through g with the tangent 
vector g are the following: 

g(e) = e eT, g, r] = R g - U g, 
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and 

9(e) = ge er > , rj = L g -i*g, 

which allows to establish the connection of the gradient V/ with the (somewhat more convenient) 
notions of the left and the right Lie derivatives of a function / : G h 1: 

Vf(g) = R* g - 1 df(g)=L* g - 1 d'f(g). 

Here df(g) :G^g* and d'f(g) :G^g* are defined via the formulas 

d 



(df(g), V ) = -f(e^g) 
(d'f{gU)=j e f{9^) 



6=0 



e=0 



Vr/ G fl , 
Vt? G g . 
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B Lagrangian equations of motion 



CONTINUOUS TIME 


DISCRETE TIME 


General Lagrangian systems 


Ho, 9) 


L(sfc,5fc+i) 


f n = VgL 
I n = v 9 L 


f n fe = -ViL(g fc ,5 fe+ i) 
1 n fe+1 = V 2 L(.g fc ,.g fe+1 ) 


Left trivialization: M = L*U 


n = Lg-i^g 


L(Sfc,5fc+i) = L (0 (5fe, W k ) 
W k = 9k V+i 


M = L*U = VqL« 


M k = L* gk U k = d{ v U l ){g k _ 1 , WVi) 


f M = ad* Tl ■ M + d'gL^ 
\ g — L 9 *fi 


f Ad* W^ 1 ■ M k+1 = M k + d'gh^ (g k ,W k ) 


Left trivialization, left symmetry reduction: M — L*H, P — Adg 1 -( 


L(g,g)=CW(P,Cl) 

n = L g - u g, P = Adg- 1 -( 


U9k,9k+i) = ^ l) (Pk,W k ) 
Wk=9 k 1 9k+i, Pk = A&g- 1 ■ C, 


M = L* g Il = Vq£W 


M k = L* g U k = d^AW(P fe _i, Wfc-i) 


f M = ad* n-M + ad* P-V P C^ 

\ P = [P,n] 


f Ad* W^ 1 ■ M k+1 = M k + ad*P k ■ V P A(')(P fc , W k ) 

\ p k+1 = kdw k - 1 -p k 
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Right trivialization: m = R*Il 


L(g,g) = L^(g,uj) 
uj = R g -i*g 

m = R* g U = V w l/ r > 

j rh = —ad* u> ■ in + dgL^ 
\ g — R g *u> 


I 
I 

r 

\ 


'(9k,9k+i) = L( r )(3 fe ,w fe ) 

Jk = 9k+l9k 

n k = R* 9k U k = d w li r \gk-\,Wk-i) 

f Ad* w k ■ m k+ i — m k + d g l^ r \g k , w k ) 
9k+i = Wk9k 


Right trivialization, right symmetry reduction: to = R*H , a — Adg ■ C 


L(g,g)=£^(a,uj) 

u = R g -i*g, a = Adg-( 

m = R g U = V U £M 

( to = —ad* to • to — ad* a ■ V 
y a = [uj, a] 


I 
i 

r 

\ 


'(9k,9k+i) = A^(a k ,w k ) 
Jk^gk+igk 1 , a k =Adg k -Q 

n k = R* k Uk = d w A^(a k -i,w k -i) 

Ad* w k ■ m k+ i =m k - ad*a fe • V a A (r) (a fe , w k ) 
L a k+1 = Adw k - a k 



The relation between the continuous time and the discrete time equations is established, if we 

set 

9k = 9, 9k+i = 9 + eg + 0{e 2 ) , h(g k , g k+1 ) = eL(g, g) + 0{e 2 ) ; 
Pk = P, W k = 1 + eU + 0(e 2 ) , A^(P k ,W k )=eC (l \P,n) + 0(e 2 ); 
a k = a, w k = 1 + euj + 0(e 2 ) , A^\a k ,w k )=eC ( - r \a,io)+0(e 2 ). 

C On 577(2) and su(2) 

The Lie group G = SU (2) consists of complex 2x2 matrices g satisfying the condition gg* = g*g = 
1, where 1 is the group unit, i.e. the 2x2 unit matrix, and * denotes the Hermitean conjugation, 
i.e. g* = g T . In components: 

_(a(3\_(a + ibc + id\ . . 

g ~ \ -0 a )-{ -c + id a-ib ) ' (LjA) 

where 

a = a + ib , [3 = c + id £ C , 
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and 

\a\ 2 + |/3| 2 = a 2 + b 2 + c 2 + d 2 = 1 . (C.2) 

The tangent space T e SU(2) is the Lie algebra g = su(2) consisting of complex 2x2 matrices rj 
such that rj + rf = 0. In components, 

V =( lb . A C+ l d ) . (C.3) 
1 y-c + id -ib J y ' 

The Lie bracket in su(2) is the usual matrix commutator. 

Let us introduce the following notations: for an arbitrary matrix g of the form ( p.l| ), not 
necessary belonging to SU(2), set 

so that $t(g) is a scalar real matrix, and ^s(g) G su(2). 

As it is always the case for matrix groups, we have for g G SU(2), ij G su(2): 

L g *r] = gr] , R g *rj = rjg , Ad g ■ rj = g^g' 1 . (C.4) 

If we write ( |C.3| ) as 

i/ = if "^."^ ) , (C.5) 

2 \ T]2- IT)1 tT] 3 J 

and put this matrix in a correspondence with the vector 

v = (m,m,V3) T gr 3 , 

then it is easy to verify that this correspondence is an isomorphism between su(2) and the Lie 
algebra (k 3 , x), 

where x stands for the vector product. This allows not to distinguish between 

vectors from R 3 and matrices from su{2). In other words, we use the following basis of the linear 
space su{2): 

1 / -i \ 1 1 / -1 \ 1 

61 = 2 ^ -z J = 27 (T1 ' 62 =2^1 J = 27" 2 ' 

-Ko 'h**' (c - 6) 

where aj are the Pauli matrices. 

We supply su(2) with the scalar product (•, •) induced from R 3 . It is easy to see that in the 
matrix form it may be represented as 

(77,C>=-2trfaC) = 2trfaC*). (C.7) 
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This scalar product allows us to identify the dual space su(2)* with su{2) itself, so that the coadjoint 
action of the algebra becomes the usual Lie bracket with minus: 

adV£ = itv] =-adr?-C. (C.8) 

We use a formula similar to ( |C.7| ) to define a scalar product of two arbitrary complex 2x2 matrices: 

( gi ,g 2 ) = 2 tr( gi g* 2 ). (C.9) 

(In particular, the square of the norm of every matrix g £ SU(2) is equal to 4). The formula (|C.S| ) 
gives us a left- and right-invariant scalar product in all tangent spaces T g G. Indeed, to see, for 
instance, the left invariance, let gr],g( £ T g SU(2) (here g € SU(2), n,( £ su(2)). Then 

{9V,g() = 2tr(OTC*5*) = 2tr(7/C*) = (77, • 

This scalar product allows to identify the cotangent spaces T*G with the tangent spaces T g G. It 
follows easily that: 

L*^ = g-^, R* g C = ^9~\ Ad* g •£ = g'^g . (CIO) 

(in these formulas g £ SU(2), so that g -1 = 5*). 

Let us now formulate several simple properties of SU(2) and su(2) which will be used later on. 



Lemma C.l For an arbitrary g £ SU(2): 

g = cos(9)l + sm(9)( with (esu(2), (£0=4. (C.ll) 

The adjoint action of SU(2) on su(2) has in these notations the following geometrical interpretation: 
grig -1 is a rotation of the vector r\ around the vector £ by the angle 29. 

This interpretation makes SU{2) very convenient for describing rotations in R 3 (in some respects 
more convenient than the standard use of SO(3) in this context). Since by rotations only the vectors 
on the rotation axis remain fixed, we see that for the case G = SU (2) 

GKl = GfW . 

In a different way, the previous lemma may be formulated as follows. 
Lemma C.2 For g £ SU{2), if 

2^f- =£, (C.12) 

then 

5 = ^M(l + e) , tr( g )= - 2 ==- (C.13) 

2 V qj ' Ky) VI + (£,0/4 

We have also the following simple connection between the matrix multiplication and the com- 
mutator in su{2) . 
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Lemma C.3 For rj,( £ su(2) their matrix product has the form ( |C.1| ), and 



(C.14) 



In particular, the following corollary is important: 



(77,0 = => v( 



Cv ■ 



(C.15) 
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